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We show that for all the three standard symmetry classes (unitary, orthogonal and symplectic), 
the conventional replica nonlinear a model gives the correct non-perturbative result for the two-level 
correlation functions Riiyi) of electrons in disordered metals in the limit of large uj. In this limit, 
non-perturbative oscillatory contributions arise from a degenerate saddle-point manifold within this 
a model which corresponds to the replica-symmetry breaking. Moreover, we demonstrate that in the 
unitary case the very same results can be extracted from the well known exact integral representation 
for Ri(iS). 

I. INTRODUCTION 

Starting from the seminal papers of WegnerEl and Efetov,i a field-theoretical description based on the nonlinear 
<7 model (NLcrM) has become one of the main analytical approaches to various problems in disordered electronic 
systems. The ensemble averaging over all configurations of disorder is performed either using bosonicalj atdermioniccl 
n- replicated fields and taking the limit in the results, or using supersymmetric (^-graded) fields.cltl 

The very first application of this approach was a derivatioruu'El of the renormalization-group (RG) equations of 
the scaling theorytj of Anderson localization. For such a perturbative derivation, generalized later for mesoscopic 
systemsEI, both the replica and the supersymmetric methods are equally well justified: they just ensure the cancellation 
of unphysical vacuum loops in a diagrammatic expansion. However, it soon became conventional wisdom that there 
existed two sets of problems for which only one of these methods was applicable. n 

On the one hand, the fermionic replica NLcrM has been generalized by FinkelsteinEl to include interactions between 
electrons. The interest in this approach has been greatly enhanced by the recent discoveryS of an apparent metaL 
insulator transition in 2D disordered systems in zero magnetic field. Although it has been recently demonstrated!!!!! 
that the Keldysh technique provides a viable alternative to the replica approach, the latter still remains one of the 
best available tools for consideration of interacting electrons in disordered systems. In the very least, it is clear that 
there is no simple way of applying the supersymmetry method to a many-particle fermionic system. 

On the other hand, the viability of the replica approach was undermined by the existence of a set of problems which 
could apparently be solved only with the. help of the supersymmetry method. The first, and possibly most famous, 
of such problems was solved by EfetovEd who used the supersymmetric NLcrM to derive the two-level correlation 
function (TLCF) in the universal ergodic regime for electrons in disordered, metallic grains. The results proved to 
be identical to those for-eigenvalue correlations in random matrix theory,E3 as had much earlier been assumed by 
Gor'kov and Eliashberg.E^I This first microscopic derivation of essentially non-perturbative results has opened the way 
to numerous new results (mainly obtained during the last decade) for which using the supersymmetry method seemed 
to be absolutely essential (many of these results have been reviewed in aj-ecent book by EfetovEJ) . 

In many ways, though, this first non-trivial "supersymmetric" resuliO seemed to be the best illustration of why 
the replica method could only be used within a perturbative approach. For the easiest case of the unitary symmetry 
the irreducible TLCF is given by 

W = -^, (id) 

where u is the distance between two levels in units of A/ it and A is the mean level spacing. This result is valid in 
the ergodic regime, i.e. for u> <C g where g 3> 1 is the dimensionless conductance. For w 3> rL, the TLCF averaged 
over fast oscillations, could be readily obtained from the standard diagrammatic techniquestJ with 1/oj being the 
perturbation parameter (or from the perturbation theory in the framework of either the supersymmetric or replica 
NLcrM). However, the non-perturbative factor sin 2 w cannot be restored from the perturbation series. Since the re plica 



trick is well justified only within the perturbative approach, it might seem rather hopeless to obtain the result (LI) 
within the replica approach. And indeed, quite involved calculations by Verbaarschot and Zirnbauerllj have shown that 
a direct application of the replica trick (using either the bosonic or fermionic NLcrM) has apparently not reproduced 



the TLCF given by Eq. (|Lj). 

Therefore, a very recent result of Kamcnev and MezarcO makes a real breakthrough in this area. Using the replica- 
symmetry breaking within an effective fermionic field theory corresponding to the large N limit of the Gaussian unitary 



ensemble of N x N random matrices, they have reproduced the non-perturbative result (1.1), albeit only in the region 



1 



u> 1. This raises hope that one might eventually apply the replica approach for obtaining non-perturbative results 
for interacting electrons in disordered systems. However, the method used in Ref. [l?] does not appear to be extendible 
to the standard NLcrM describing electrons in a random potential. Although the large- N held theory written in 
terms of the n x n replica matrix fields Q is, in principle, totally analogous to the standard NLcrM, this analogy has 
been considerably smeared by the very essence of the method applied. Namely, the standard exp[iVtr ln(e — iQ)} 
term has been represented after a shift of variables as (detQ)^ and the resulting Itzykson-Zuber integral (similar 
to that introduced in the context of the supersymmetric method by GuhiiiS) was calculated taking into account 
replica-symmetry breaking.EZI 

Although the method used in Ref. [l7] could hardly be generalized directly, the message was very clear: the replica 
approach taken together with the symmetry-breaking, can reproduce the results of the supersymmetric approach. 

The purpose of this paper is to demonstrate that non-perturbative results for all the classical symmetry ensembles 
(unitary, orthogonal and symplectic) can be derived in the framework of the standard fermionic NLcrM with replica- 
symmetry breaking. We believe that the present derivation is both more general and considerably easier than the 
original oneO and can be straightforwardly extended to include interactions. Independently, Alex Kamenev and Marc 
Mezard have presented in their second papeilij a derivation which is similar in spirit, although considerably different 
from that suggested here. 

Before describing our derivation, we want to emphasize the following. It is the standard fermionic NLcrM which 
does contain correct (at least for large w) non-perturbative oscillatory contributions to the level correlation functions. 
The replica-symmetry breaking which involves a set of additional saddle-point submanifolds within the standard 
NLcrM is just a very convenient method of calculating the large-w asymptotic behavior of the correlation ftmctions 
including these oscillatory contributions. Similarly, within the supersymmetric NLcrM Andreev and AltshulerEll have 
reproduced the large-w limit of Eq. (1.1) by the saddle-point calculation of the supersymmetric integral that involved 
an additional saddle-point which breaks the supersymmetry. The very same integral has been exactly calculated by 
EfetovO without breaking the supersymmetry. We extend this analogy to the fermionic replica NLcrM by showing, 
with the help of a trick which resembles the replica-symmetry breaking, that the n = limit of the exact integral 
representation of the TLCF obtained by Verbaarschot and ZirnbaueiO does contain, at least for large u, the correct 
oscillatory behavior. 

The paper is organised as follows. In the next section we recall the formulation of the standard replica fermionic 
NLcrM and introduce all relevant notations. In Section [II we present in detail a calculation of the large ui non- 
perturbative contribution to the TLCF for the unitary ensemble, based on replica-symmetry breaking. In Section 
IV we outline similar calculations for the orthogonal and symplectic ensembles. In Section ^ wft_extract the very 
same results from the exact integral representation for R 2 obtained by Verbaaschot and ZirnbauerEj for the unitary 
ensemble. Finally, in Section VI we summarize our results and discuss perspectives of the replica method. 



II. FORMULATION OF THE PROBLEM 



We consider the two-level correlation function (TLCF) defined by 



-R2M = - s (v(e + u)v(e))-l. 



(2.1) 



Here (...) stands for the ensemble averaging, v(e) is the electronic density of states per unit volume defined in terms 
of the spectrum {e a } for a given sample as v{e) = L~ d J2 a S( £ ~ £ a), v = — l/£ d A, and L, d are the sample 

size and dimensionality. We will measure all energies in units of A/-7T, in which the TLCF is as written in Eq. flLl| ). 
In these units we express R 2 via the product of the retarded and advanced Green's functions as follows: 



The function S 2 can be expressed in the standard way0'0 in terms of the replica NLcrM: 

1 d 2 



lim 



n—>o n 2 duj 2 

where T[Q]oS\ is the non- linear sigma model functional 

1 



Z n (w), Z n (u)= VQexp(-F[QM) 



T\QM = j- d I dVrr 



\D(y Q ) 2 - ^AQ 



(2.2) 



(2.3) 



(2.4) 



2 



This functional has been derived^ for the model of free electrons in a random potential represented in terms of the 
fermionic (anticommuting) replica fields. As a result, the field Q(r) is an Hermitian matrix of rank 2n whose elements 
are quarternions for the orthogonal symmetry ensemble, complex numbers for the unitary symmetry ensemble which 
arises when time-reversal symmetry in the model is broken, and real numbers for the symplectic symmetry ensemble, 
which arises when spin-rotation symmetry is broken. It satisfies the standard constraints 

Q 2 = t 2n , Tr Q = . (2.5) 

These constraints are resolved by representing Q as follows: 

Q = U^AU = TUT , A = diag(l n , -1„) , T = exp ( °? t * ) . (2.6) 



v -ft 0„ / 

Here t is an arbitrary n x n matrix, and T is obtained by factorizing matrices U G S (2n) with respect to redundant 
matrices R £ S (n) x S (n) that commute with A, i.e. U = RT, where S (n.)js the appropriate symmetry group. These 
conditions mean that Q belongs to the compact Grassmannian manifoldEa (coset space), S(2n)/S(n) x S(n). For 
three different classes of symmetry, S (n) becomes the symplectic group, Sp(ro), for the orthogonal class, the unitary 
group U (n) for the unitary class, and the orthogonal group, SO (n), for the symplectic class, 

The ergodic regime corresponds to the level separations much smaller than the Thouless energy, D / L 2 (which in the 
chosen units coin cides, up to a numerical factor, with the dimensionless conductance g). In this regirae the gradient 
term in Eq. (2.4) may be neglected, and the NLcrM functional reduces to the zero-dimensional limitO 



T[Q-u] = - l -^Tr [AQ] , (2.7) 



with Q becoming a spatially homogeneous matrix. Here and in Eq. (2.4) a = 1 for the orthogonal class, and a = 2 for 
the unitarjj-.and symplectic classes. This factor arises because unitary and symplectic classes have been obtained from 
orthogonalo by the suppression of massive modes corresponding to the time-reversal or spin-rotational symmetry 
breaking and a subsequent reduction of the Q matrix rank. The coefficient a also absorbs an extra factor in the 
symplectic case due to the redefinition of the mean level spacing A in the chosen units. 



III. UNITARY ENSEMBLE 



We limit all further considerations to the ergodic regime only, and rewrite explicitly Z(u>) given by Eqs. (2.3) and 



(2-7) 



Z n (ui) 



T>Q exp 



where the measure is defined by 



-i—TrAQ 



(3.1) 



VQ 



ra inra * 



dn = dT- t~ 



(3.2) 



Here T is the matrix parameterizing Q, Eq. (2.6), and r and a refer to the replica indices which originate from G r 
and G a , respectively. In the large-w limit this integral is mainly contributed by the extrema of the functional which 
obey the standard condition [A, Q] = 0. This condition is satisfied by any matrix of the form Q = diag(Q r , Q a ), 
where Q r and Q a are the n x n Hermitian matrices whose eigenvalues are ±1 and Tr (Q r + Q a ) = 0. This defines a 
highly degenerate saddle-point manifold which consists of submanifolds specified by a particular distribution of 
n eigenvalues '+1' and n eigenvalues '-1' between Q r and Q a . These submanifolds can be divided into n+ 1 classes of 
equivalence, Q p — diag(Qp, Q^), labeled by TrQ^ = — TrQ° = n — 2p, with p = 0, 1, . . .n. The p-th class has weight 
(C£) 2 ,withC^Q. 

The matrix with (Qp) 2 = l n and TrQ^ = n—2p can be parameterized by analogy with Eq. ( p"7j| ) as 



Q r p = (T;?X p T; , A p = diag(l„_ p , -l p ) , T; = exp 



t r 

-(r)t o P 



(3.3) 



where t r is an arbitrary p x (n—p) matrix. This defines the coset space G p = U {n)/\J (n—p)x U (p), i.e. Therefore, 
Q p = diag(Qp, Q°) belongs to the manifold G p x G p and can be parameterized as 
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Qp T p A p T p ! 



Tp = dia g (r;,T;) 



A. 



diag(Ap, -Ap) 



(3.4) 



The integer p specifies the replica-symmetry breaking, as it describes the number of the —1 eigenvalues in each Q r 
block (equal to the number of the +1 eigenvalues in each Q a block): in the symmetry- unbroken case, p = 0, and 
hence retarded and advanced blocks, Q r,a , contain only positive or negative eigenvalues, respectively. We want to 
emphasize that m anifo lds G p x G p ' with p =^ p' cannot appear within the degenerate saddle-point manifold, [A, Q] = 0, 
of the functional (2.7). Indeed, the corresponding matrices have non-zero trace, and thus do not belong to the Q 
space defined by Eq. (2.6). Naturally, one could have derived the replica NLcrM with different numbers of elements 



corresponding to the retarded and advanced blocks, as in Ref . |. In this case, the replica-symmetry breaking would still 
mean the redistribution of +l's and — l's eigenvalues between these blocks, keeping the value of Tr Q the same (equal 
to the difference between the numbers or replicas in the retarded and advanced blocks in the the symmetry-unbroken 
case, p — 0) in each p-th manifold. 

Now we need to take into account contributions from 'massive' modes (with mass oc 1/lu) in the vicinity of each 
manifold (3.4). In the large-w limit these contributions may be considered as independent and the partition function 



is then represented by the sum of all of them: 



Z n {uj) = 



£(C£) 
p=0 



PQexp — i^TrA 



(3.5) 



Here we have used the fact that the similarity transformation in the vicinity of A p , i.e. UA p U^, covers the entire 



symmetric manifold of the NLcrM, Eq. (2.6), including all the massive modes. Having substituted this into the 
functional fl2.7p , we have reduced Tr[/A p [/ T A to Tr A p T^AT = TrA p Q, where we have substituted U = RT, as 
defined after Eq. (2.6). Let us emphasize again that the above representation of Z as the sum over all p can be 
justified only as a perturbative (in l/u>) procedure: a possible overlapping of massive modes originated from different 
manifolds is irrelevant in the large-w limit. 

Each term in the sum ( |3.5| ) contains both massive and massless modes. Indeed, we have used above the factorization 
U = RT with R being block-diagonal matrices commuting with A. The matrices T in TrA p Q = TrATApT''' still 
contain the subset of matrices commuting with A p that correspond to the massless modes. Therefore, we need to 
parameterize T in a way which enables us to factorize out these massless modes and perform the integration over the 
massive ones. 

The most suitable parameterization of T, analogous to that used in Ref. [l6] for the bosonic NLcrM, can be obtained 
by expanding the matrix exponent in Eq. ( [2.6| ). By introducing matrix B = t(tH)~ 1 ^ 2 smVtH, we represent T and 
thus Q = TUT as follows: 



T 



y/l n - BB^ B 

-B 1 y/t n - B^B 



Q 



- 2BB^ B^t n - BiB \ 
BV*n - B& -(l„-2St J B) J 



1, 



(3.6) 



The matrix B in this parameterization is not unconstrained, though. The Q — condition is fulfilled only when 
the matrices \/l n — BB^ and y/l n — B^B are Hermitian. This is so only when a ll th e eigenvalues of BB^ and B^B 
do not exceed unity. Only under this constraint does Q, parameterized as in Eq.SJ), still belong to the coset space 
U(2n)/U(n) x U (n). Nevertheless, this parameterizatio n is very convenient. First, the corresponding Jacobian is 
unity (see Appendix), so that the measure of integration (3.2) can be written simply as 



(3.7) 



In addition, the representation of all the exponents in the sum (3.5) in terms of B is also very simple, TrA p Q 
2(n - 2p) - 2Tr \ p (BB^ + B^B), so that we obtain: 



Z n H = £(C£) 2 -e^-») Z P {uh 



p=0 



ZP(uj) = J VBcxp [iuoTrXpiBB^ + B*B)] 



(3.8) 
(3.9) 



The region of integration in (3.9) is restricted by the constraint described after Eq. ( |3.6| ). Last, but not least, the 
parameterization (3.6) allows us to separate out the massless modes, which obey the condition [T, Ap] = 0, in each 



4 



integral (3.9). Indeed, this condition is satisfied by all matrices T constructed from B which anticommute with A p , 
i.e. have the off-diagonal block structure. 

This means that in the representation of B in the block form reflecting the structure of \ p — diag(l^ 



L n—p: lp)j 



B = 



Bx h 
4 B 2 



(3.10) 



the matrices B\ 2 represent the massiv e m odes, and b\ 2 massless. When the massive modes are suppressed (B\ = 
and B 2 = 0), the T matrices in Eq. (3.6) constructed from p x (n— p) matrices 61.2 only, parameterize the same 
degenerate p-the manifold, G p x G p de scribe d in Eq. ( |3.4|) , as one expects. 

By substituting the representation ( 3.10 ) into Eq. (|3.9|), we reduce Z% to the product of integrals over the massive 
and massless modes: 



Z*(w) = / VB X VB 2 exp f-2wjtr {B x B\ - B 2 B\ 



T>biVb 2 , 



(3.11) 



Here the region of integration over 61.2 d epends on B\ 2 due to the constraint on the eigenvalues of the matrices 
BB^ and B^B in the representation ( 3.10| ). Since the integral over B\$ is contributed only by the region where both 
tr B\B\ and tr B 2 B\ < 1/lu -C 1, in the leading in 1/ui approximation we may put both B\ 2 to in the cons train t 
of the integration region over the massless m ode s &i,2- In this approximation, as we have noticed after Eq. (3.10), 
matrices 61.2 parameterize the p-ih manifold (3.4) so that 



Jvb 1 Vb 2 = JvQ p = n 2 (G p ) 



(3.12) 



where the measure of integration over T>Q p is defined in terms of T p in the same way that T>Q is defined in terms of 
T, Eq. (3.2), and f2(G p ) is the volume of the compact coset space G p . This volume is expressed via the well-known 
volumes of the unitary group, f2(U (n)), as follows: 



n(G p ) = 



fi(U(n)) 



fi(U(n-p))0(U(p)) 



= (2tt) 



L[« 2 -(«- P ) 2 ^ 2 ] 



n 



r(i + j) 

lr(n + 2-j)' 



(3.13) 



In the same large-w approximation, the variables B\ 2 parameterizing the massive modes are unconstrained. Then 
the Gaussian integral over the 2[(n—p) 2 + p 2 ] independent massive modes yields 



VBx VB 2 exp 



2iuj\x{B 1 B\ - B 2 B\ 



\10J/ 



(3.14) 



Combining Eqs. ( 3.14 ) and ( 3.13 ) and omitting an irrelevant overall factor which goes to 1 when n — > 0, we arrive at 
the following expression which is essentially the same as that derived in Ref. [l^ via the Itzykson-Zuber integral: 



p=0 



e «j(2p-n) 

(2cj)(™ _ p) 2 +p 2 



k = ci n 



r(i + i) 



JLi r( n + 2 - j) 



(3.15) 



Here the summation over p has been extended t o 00 s ince F% — for all integer n > p. This allows one to take the 
replica limit, n — > 0, in each of the terms in Eq. ( 3.15 ) . D ue to the fact that fv « p 



p = and p = 1 in Eq. (3.15) contribute to S 2 in Eq. (2.3). Omitting all the terms with p > 2, one obtains 

g — tun g2uj(2 — n) 

= — WT- 



4cj("- 1 ) 2 + 1 



oc rt F as n — > 0, only the terms with 

(3.16) 



Substituting this into Eqs. ( |2.3| ) and ( |2.2| ) and keeping the leading in lw terms only, one arrives at the expression (1.1) 
for the TLCF. Although this expression looks as being exact, it has been actually derived only in the large-w limit, as 
is the case in the Itzykson-iZuber calculation of Kamenev and MezardE3, and in the 'supersymmetry breaking' method 
of Andreev and Altshulcr.E] 
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IV. ORTHOGONAL AND SYMPLECTIC ENSEMBLES 



A generalization to both these symmetries from the unitary one is straightforward' .elements of all matrices in the 
previous section were complex numbers, while now they become quarternion numberstjt 2 ] in the orthogonal symmetry 
and real numbers in the symplectic sym metry. In other words, S(n), S(p), etc, should be substituted for U (n), 
U (p), etc, in all formulae of Section III, with S being the symplectic group Sp for the orthogonal symmetry and 
the orthogonal group SO for the symplectic symmetry (such an inve rsio n is characteristic of the fermionic replica 
approach). One should also redefine the m easu re of integration, Eq. (3.2), in terms of independent matrix elements 
and substitute 2uj by acu as defined in Eq. (2.7). 

After these redefinitions, the partition function Z n is represented by the sum ( |3.g| ), with each term in this sum 
given by the product of integrals over the massive and massless modes, Eq. ( 3.11 ). The number of massive modes is 
(4//3)[(n — p) 2 + p 2 ], where j3 = 1, 24 for the orthogonal, unitary and sympl ectic symmetry is the standard Dyson 
parameter. For orthogonal and symplectic symmetries, the Gaussian integral (3.14) over the massive modes yields 



(tt/lu) 21 ^ 2 ^ (-1)(«^) 2 -p 2 orthogonal symmetry 
(vr/2 W ) [( "^ )2+p2l/2 i -[(«^) 2 -P 2 l/2 symplectic symmetry 



(4.1) 



The integrals over massless modes, Eq. ( 3.12 ), give the volumes, fi 2 (Gp) of the corresponding coset spaces. The results 
can be written as 



r(i + 2j) 



r(i + j/2) 



(7T 



(n-p)p/2 



n 



Lr[l + (n+l-j)/2] 



orthogonal symmetry 
orthogonal symmetry 



(4.2) 



Combining Eqs. 4.1 and 4.2, we obtain the partition function for both the symmetries. It can be written in terms of 



/?, together with the unitary Z n , Eq. (3.15), in universal form: 

— iu{2p— n) 



p=0 



[(n-p) 2 +p 2 



r(i + fj) 



=i r (! + f (^ + 1 -.7)) 



(4.3) 



In the n — > limit only the terms with p = 0, 1 survive for /3 = 1, 2, while for f3 = 4 one needs to keep the p = 2 term 
as well: 



Z n ^{uj) = 1 — iujn + n 2 < — 



r 2 (i + 2//3) 



Substituting this into Eqs. 



and (2.2), we obtain the TLCF as follows: 
1 



Ra(u) 



f3u) 2 



2T 2 (l + 2/f3) n . cos4w 

COS 2UJ + O0A- 



32uj 4 



(4.4) 



(4.5) 



This reproduces the correct large-w asymptotic behavior of R2 for all three ensembles, which includes the non- 
perturbative oscillatory factors. 

As we have derived these results with a non-standard treatment of the standard fermionic NLtrM, they also should 
be contained in the exact integral representation obtained for this model by Verbaarschot and ZirnbauerEj We will 
show now that this is, indeed, the case. 



V. THE LARGE-W LIMIT OF THE VERB A ARSCHOT-ZIRNB ALTER INTEGRAL 



The 'zero-mode' partition function, Eq. (2.7), can be exactly represented in the following form 
Z n (oj) = n 2 ((J (n)) f +1 A 2 (A) f[ e-^dXi , A(A) = - A,) . 

i=l i<7 



(5.1) 



G 



which is equivalent (with accuracy up to factors going to. 1 in the n — > limit) to the representation for S2 given in 
Eq. (2.24) of the paper by Verbaarschot and Zirnbauerilj 

The leading in contributions to this highly oscillatory integral (which does not have stationary points inside 
integration region) come from the end points. To single out these contributions, we must take some A's close to +1 
and the rest close to — 1, which imitates replica symmetry breaking. Let us choose n—p of A's close to +1 and p of 
A's close to —1. Then we can split up the Vandermonde determinant in the following way: 

n 

A 2 (A) = Y[ \Xj - Ai| pa 2 2p( - n ^ A 2 + A 2 _ 

where 

n~p n 

A 2 + = J] |A«-A,-| A 2 _ = [] ( 5 - 2 ) 



Reducing the integral (5.1) to the sum of such contributions only, we represent it as 

n . + 1 n-p +00 n 

Z n (uj)^n 2 ((J(n))^2{CP)2 2p ( n - p '> A 2 + Y[e- iuX * dAj / A 2 _ J] e"^ dX, . (5.3) 

p=l j=l J ~ X j=n-p+l 

Since in each of the integrals all the variables are close to one of the limits of integration, the second limit was extended 
to infinity. Now we make substitutions Aj = 1 — Xi in the first integral, and A^ = — 1 + SEj in the second one, reducing 
the above sum to the form: 

n 

Z n {uj) » fi 2 (U (n)) £(C*) 2 2p ^- p h^ 2 ^I n ^I p (5.4) 
P =i 

where I p are integrals of the Selberg's type:@ 

Substituting the known Selberg integrals and discarding an overall factor which goes to unity in the replica limit we 
arrive at 

Z n (co) = V [F p ] 2 ■ —. w -^, F'P = CI TT w 1 Z ^ ^ (5-6) 

p=0 j'=l v Jy 



This expression is exactly the same as Eq. p. 15 obtained in Section HI with the help of the replica-symmetry breaking 



Therefore, the exact representation (5.1) does contain the true oscillatory asymptotic behavior of the TLCF. 

The authors of Ref. |l6| have also drawn attention to the fact that there is an apparent contradiction between the 
lu = limit for S2 obtained from the replica trick and the exact supersymmetric result. Indeed, if ui is put to in the 



expression for following from Eqs. (5.1) and (2.3) the n — > limit is taken after that, one obtains S 2(^ = 0) = — 1. 
This cannot be correct as 3?e5a(w — > 0) > as follows from the definition ( |2.2| ) and, moreover, it is known that 
£2(0;— *0) — * 6 (u>). What is interesting, however, is that if one separates the singular, ^"^(w), an d regular, S™ s (u)), 
parts of S , 2(w), then Sl es (uj—>0) = —1. Therefore, the replica method gives S' 2 cs (u; = 0) correctly, and it is just S^ 1 " 8 
which is missing. However, the fact that Z n (uj — > 0) is finite for any integer n does not necessarily implies that it 
is also finite (as a function of u> — > 0) in the replica limit. For example, if the expansion of S^w) before taking the 
n = limit contained a term proportional to ui n , it would be singular in the replica limit. In other words, if a 
non-trivial dependence on the order of limits n — > and w->0 existed, in a spirit of the replica trick the n — > limit 
should be ta ken first. At the moment, though, this is only a speculation as we have not yet succeeded in calculating 
t he integral (fO) for small to. However, the fact that the large-w limit of this integral reproduces the correct results 



(1.1) makes it plausible that there are only technical difficulties rather than one of principle in the application of the 



replica method. 
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VI. CONCLUSIONS 



The purpose of this paper was to demonstrate explicitely that non-perturbative oscillatory contributions to the 
TLCF of electrons in a random potential could be extracted from the standard NLcrM formulated in fermionic 
replicas by Efetov, Larkin and Khmclnitskiicl many years ago. To this end, all one needs is to parameterize all the 
non-trivial saddle-point manifolds corresponding to the broken replica-symmetry and describing 'massless modes' of 
the theory, and expand the action in the vicinity of these manifolds to include 'massive modes'. The very sumilar 
approach has been used in the supersymmetric NLcrM: the non-perturbative oscillations have been extracted by the 
expansion around two extremal points one, of which breaks the supersymmetry.EEl Since the, exact supersymmetric 
calculation of the TLCF was well knownjli it was clear that the supersymmetry breakingpl was just a convenient 
method of extracting the large-w limit (and going beyond the universal 'zero- mode' apporximation) . It is not clear 
at the moment whether exact calculation is possible within the replica approach. However, as we have shown in the 
Section |v|, the exact integral representation of i?2(w) does at least describe the correct behavior in the large-w limit. 

On the other hand, there is a number of non-perturbative results where all the leading contributions come from the 
small uj limit (see, e.g., Ref. |23| ): the existence of a nontrivial spatially inhomogeneous saddle point in the non-compact 
sector of the supersymmetric NLcrM was the main source of these contributions. At the moment it is unclear whether 
such a saddle point could be found-within the fermionic replica approach. Nevertheless, the fact that some of these 
results had earlier been obtainedLrO within the RG treatment of the fermionic replica NLcrM prompts that it could 
be technical rather than methodoligical difficulties involved. 

Despite the fact that the replicas prove now to be a trusted tool for non-perturbative calculations, there is no doubt 
that the supersymmetric NLcrM is the best tool to obtain non-perturbative results for non-interacting electrons. But 
the further elaboration of the non-perturbative methods within the replica trick will be important for the problems 
where supersymmetry cannot be applied. This is the situation with the interacting electrons in a random potential. 
Such systems can be described by NLcrM with fermionic replicas only and we believe that the further development of 
methods similar to presented here will enable one to obtain at last non-perturbative results for disordered systems of 
interacting electrons. 
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APPENDIX A: CALCULATION OF JACOBIAN 



To calculate differentials of matrices entering the measure (|j) we use the following rational substitution: 

1 + rrt 2 

2 / 

B 



1 - rrt 



T 



1 — tt^ 1 — tt^ 



1 -tV 



4 



1 + T T T 
1 - T+T 



Then the differentials in the measure ( |3.2| ) are given by 

dft 12 (T) = (dT • T f ) = 2 — r (dr -rdrr) L- 

1 — TT T 1 — T T 7 



On the other hand we can find the differential of B, entering the measure (3.7) as follows: 

1 .si 



dB 



!- r (dr + r 

1 — T"H 



dr r) 



1 - tN 



(A.l) 



Comparing Eqs. ( A. 2 ) and (A.3) one can see that the Jacobians of transitions dfl 12 (T) — 
exactly the same. From that fact one deduces that the Jacobian of the transition d£l 12 (T) 



dr and dB — 
-> dB is unity. 



(A.2) 

(A.3) 
dr are 
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